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We describe the statistics of temperature fluctuations in a SINIS structure, where a normal 
metal island (N) is coupled by tunnel junctions (I) to two superconducting leads (S). We specify 
conditions under which this structure exhibits manifestly non-Gaussian fluctuations of temperature. 
We consider both the Gaussian and non-Gaussian regimes of these fluctuations, and the current 
fluctuations that are caused by the fluctuating temperature. We also describe a measurement setup 
that could be used to observe the temperature fluctuations. 
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I. INTRODUCTION 

Nanoelectronic structures with superconductor- 
insulator-normal metal (SIN) junctions are often used in 
applications involving thermometry or cooling,^ enabled 
by the existence of an energy gap of 2A in the density of 
states of the superconductor. In those systems, typically 
two such junctions are connected in series to form a 
SINIS structure. Thermometry with such structures 
relies on the temperature dependence of the electric 
current through the SIN junction,^ whereas refrigeration 
is based on the heat current carried by electrons from 
the normal metallic part.'^''* 

In most SINIS devices the electron-electron relaxation 
rate is fast compared to the electron injection rate due 
to driving, and therefore the electron system in the de- 
vice maintains its equilibrium energy distribution, de- 
scribed by the Fermi function. In this quasiequilibrium 
case the temperature of the electron system may fluctu- 
ate around its average value:^ The average temperature is 
determined by a heat balance, where incoming and outgo- 
ing heat flows cancel each other. These heat flows exhibit 
fluctuations due to the probabilistic nature of the tran- 
mission of electrons through the tunnel barriers formed 
by the insulating layers. As a result, the temperature of 
the electron system fluctuates as well. 

Recently, temperature fluctuation statistics (TFS) has 
been studied in non-interacting islands,^ and overheated 
single-electron transistors.®'^ The TFS become especially 
interesting in SINIS structures, where a normal metal- 
lic island with a single-particle level spacing 5 ^ A is 
connected to two superconducting reservoirs via tunnel 
junctions.® By virtue of the interplay of regular quasipar- 
ticle tunneling and two-electron Andreev tunneling, the 
statistics of temperature fluctuations are strongly non- 
Gaussian when all transmission eigenvalues of the junc- 
tions are small, characterized by a number < ((5/A)^/^, 
and when the voltage is close to 2 A, \V — 2A| < ^/5~K. 

In this Article we apply the theoretical methods devel- 
oped in Ref. 7 to the SINIS structure, shown schemati- 
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FIG. 1. (color online) Schematic diagram of a SINIS struc- 
ture, biased with voltage V . The normal metallic island is 
connected via tunnel contacts to superconducting reservoirs 
at temperature Ts- The transmission eigenvalue distribution 
of the contacts is characterized by two numbers, g and J-. The 
temperature on the island, T, fluctuates due to fluctuations 
in the energy flows to the leads and the phonons. Tunnel- 
ing through the contacts is described by the connector action 
Sc, and the electron-phonon coupling with the corresponding 
action 5e_ph. 



cally in Fig. 1. We present a detailed theoretical descrip- 
tion of the temperature fluctuations in the Gaussian and 
non-Gaussian regimes and show how these fluctuations 
give rise to large fluctuations in the electric current, with 
a noise power that can exceed the intrinsic current fluctu- 
ations by a factor of a hundred. We extend the previous 
theoretical description to include electron-phonon cou- 
pling and show that its effect on our results is negligible, 
provided that the superconducting reservoirs can be kept 
at a relatively low temperature, Tg < O.IA. 

This Article is structured as follows: Section II outlines 
the theoretical methods utilized for the calculation of the 
TFS. In Sec. Ill we focus on the temperature and current 
fluctuations in the Gaussian regime and in Sec. IV on the 
fluctuations in the non-Gaussian regime. We describe a 
method for the detection of these fluctuations in Sec. V. 
We conclude and discuss our findings in Sec. VI. In the 
entirety of the text, we use units in which h = ks = e = 
1. 
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II. EFFECTIVE ACTION FORMULATION 

We study the temperature fluctuations using the 
Keldysh action technique, reviewed in Ref. 9 and de- 
tailed for the application to temperature fluctuations in 
Refs. 5 and 7. We assume that the reservoirs are diffusive 
bulk superconductors in equilibrium, and that the nor- 
mal metal island is also diffusive and in quasiequilibrium. 
The condition for the island to be in quasiequilibrium is 
given by {G /GQ){Gi/GQ){6 / ma^{T,V)} < 1, where G 
is the total conductance through the structure, Gi is the 
conductance of the island not including the tunnel barri- 
ers, and 5 is the single-particle level spacing on the island. 
The conductance quantum is defined as Gq = 2e^ /h. 

After incorporating the conservation of energy E and 
charge Q on the normal metal island with the aid of two 
Lagrange multipliers, ^ and x/*^ the Keldysh partition 
function for the SINIS structure can be written in the 
form 

Z= j V^{t)Vx{t)VE{t)VQ{t)expl^- j dt {^{t)E{t) 
+x{t)Q{t)-S{i,X,E))], (1) 



where S = Sc + >Se-ph is the action of the system, con- 
sisting of the connector action Sc, describing the con- 
tacts to the superconducting reservoirs, and the electron- 
phonon action (Se-phi describing the coupling of electrons 
on the island to the lattice phonons. For the total energy 
and charge on the island we use a simple model where 
E = 7r2TV(6(5) and Q = C^i. Here, T is the tempera- 
ture, ji the chemical potential, and C the capacitance of 
the island. The single-particle level spacing depends on 
the size of the island V, and the density of states at the 
Fermi level via 5 = l/(i^V). For a small copper island, 
V « 10^'' J~^m""^, V w 10^^^ m'^, and aluminum leads, 
A « 10-23 J, resulting in (5/A « lO'^. 

An equation describing the time evolution of the prob- 
ability distribution for E can be obtained from Eq. (1) 
in the form 



dP{E,t) 
dt 



HV{E,t), 



(2) 



where the operator H is obtained from S{^,E) with the 
substitution ^ —d/dE. The operator H must be nor- 
mally ordered, i.e., ^ must be to the left of all EJ The 
stationary solution of this equation yields the statistics 
of the fluctuations of E, and is therefore our main object 
of interest. 



A. Connector action 

The connector action for a SINIS hybrid structure is 

very similar to that of an all-normal structure.^ The 
gauge transformed quasiclassical Green's functions in 



Keldysh (") (S) Nambu (') space are given by 
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where are Pauli matrices, h = diag(l — 2/(e), 2/(— e) — 
1), and /(e) is the Fermi function. 

In terms of these functions the connector action is 
given by^^ 



de 
2^ 



Trln 



,{Ga^Gi} — 2 



(8) 



where Tr denotes a trace over Keldysh and Nambu in- 
dices, and T" are the transmission eigenvalues of contact 
a ^ L,R. When some parts of the structure are super- 
conducting, this quantity becomes relatively complicated 
due to the matrix structure in Nambu space. In the tun- 
neling limit T" ^ 1, however, this expression can be 
considerably simplified. Expanded to the second order 
in the transmission eigenvalues it reads 



5c=^E^«/^Tr[{G„,G,}-2] 

a 

-^E5«-^<^/^Tr[({G„,G,}-2 



(9) 



where ga = X^n -^n ^ 1 is the dimensionless conduc- 
tance of contact a and J"a = 'Eni^n f/T.nTn < 1 
a parameter describing the distribution of transmission 
eigenvalues. Typical value is — lO-^- lQ-^ for alu- 
minum junctions. ^2 xhe expansion to the second order is 
necessary since the first order term is exponentially sup- 
pressed at bias voltages below V = 2A. The first term 
describes tunneling of quasiparticles through the insulat- 
ing barrier, whereas the second term describes tunnel- 
ing of Cooper pairs (Andreev tunneling) and processes 
where two quasiparticles are transmitted simultaneously 
through the barrier. 

In the incoherent limit, neglecting the proximity effect 
on the island, the connector action for the left and right 
contacts can be evaluated independently, in such a way 
that the superconductor always has zero chemical poten- 
tial and the island has /it = —Hl,r for the left and right 
contact, respectively, and V = I-J.^ — iifj. This circum- 
vents the need to take into account the rotating phase of 
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the superconducting order parameter in superconductors 
at a finite voltage. 

We proceed by calculating the connector action in two 
parts. Carrying out the trace in the first term in Eq. (9) 
results in 

a 

X [2(e'^+?(^+''.)_i)/^(,)(i_/,(, + ^^)) 

+2(e-^-«(^+^=) - l)fAe + M„)(l - Ue)) . 

(10) 



Here, Re [5(e)] = Re[\e\/Ve^ - is the BCS den- 
sity of states. The second term in Eq. (9) contains 
four kinds of terms: Terms proportional to Re [(7(e)], 
(Re[ff(£)]) ^ (Re[f (£)])^ and |i^(e)|2, where F(e) = 
Asgn(e)/-\/e2 — A^. The first three are non-zero only at 
energies |e| > A, and can be neglected in comparison to 
Eq. (10) due to the small prefactor Ta, whereas terms of 
the last kind are non-zero everywhere and therefore im- 
portant. Disregarding the other terms brings the second 
term in Eq. (9) into the form 
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FIG. 2. (color online) Ratio of the voltage drops across left 
and right junctions as a function of conductance and J- asym- 
metries for V — 2A. Note the logarithmic scales. 



and are the energy current due to quasiparticle tun- 
neling through contact a and its noise. For zero temper- 
ature leads, terms in si^^ can be evaluated analytically 
with the aid of formulas from Ref. 13. For example, when 
(A + V/2)/T > (A - V/2)/T, 



-f(e-2x-2«Mo 



X [(e- 



l)(l-//(e + /i„))(l -//(-£ + Ma)) 

l)//(e + Ma)//(-e + Ma)] . (11) 



Using the connector action we can evaluate the charge 
and heat current, / = dxSc\x=o, H = d^Sc\^=o, and their 
noise powers Si = 9^5c|x=0; Si 
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For simplicity, we consider a symmetric structure with 
5l = 5fl = 5 and J-^ — J-r = T ^ and therefore 
/ii = —[iR = V/2. This is not a restrictive assumption: 
Moderate asymmetries in the properties of the junctions 
does not affect our results near V = 2 A. The reason 
for this is the strong nonlinearity of the current-voltage 
relation in the vicinity of V = 2 A, which guarantees a 
small asymmetry in the voltage drops across the junc- 
tions (see Fig. 2). We also assume that the electric RC- 
time of the structure is small compared to the character- 
istic time scale for energy transport. In this case the evo- 
lution of x{t) and ii{t) are well-described by their saddle 
point solutions,^ and when the structure is symmetric, 
x{t) ~ 0. Furthermore, we can neglect cumulants of the 
third order and higher in the quasiparticle heat current, 
since they are small at low temperatures. For this rea- 
son we expand the first part of the connector action to 
the second order in ^. The total connector action be- 



comes Sr 



5, 



(2) 



where Ha^ 



where gi{a) = — Li;(— e~°), Li; (a;) being the Zth order 
polylogarithm, and a — {A — V/2)/T. For leads at finite 
temperature T5, the heat current includes an additional 
term 



2g j 2nTs 



exp(-A/rs)A2 



(14) 



Due to the exponential dependence on Tg we neglect this 
term in the following, and thus our zero-temperature ap- 
proximation is valid as long as Ts ^ A. 

At zero temperature Sc can be evaluated to 



Sf^ =5^"— arctanh(\//2A) 
27r 

~ qJ^ — In — — 

^ 47r [l-y/(2A) 



1 



(e«^-l). (15) 



Finite temperature gives rise to logarithmic corrections 
which can be disregarded in the following. As seen from 
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the second line, the action has a logarithmic divergence 
at the threshold V = 2A, an artifact of the expansion of 
the action in transmission eigenvalues. A careful analysis 
shows that this divergence is cut off in such a way that 
the argument of the logarithm goes to a,t V = 2 A. 
However, the exact functional form of the action near the 
threshold is not necessary. Instead, the strength of the 
Andreev tunneling can be described by a single parame- 
ter in the form = Ta (e«^ - l). 



III. GAUSSIAN REGIME 

The temperature fluctuations are Gaussian when S <C 
T'^^^A or — 2A| ^ V 6 A. In this case the temperature 
distribution is characterized simply by its average value 
and variance, which can be solved from the heat current 
and its noise in a straightforward manner. 

A. Average temperature 



C. Electron— phonon action 

Elecron-phonon coupling in the lowest order is de- 
scribed with the action (see Appendix A) 



-^e— ph — 



/ dee=^[(e«^-l)n(e,Tph)(l + n(e,T)) 

(16) 



24C(5) 

+ (e-«^-l) n(e,T)(l + n(e,Tph))] 



where n(e, T) is the Bosc distribution function, S the 
electron phonon coupling constant, V the volume of the 
island, and C(5) » 1-04. The electron-phonon coupling 
strength can be adjusted by a dimensionless parameter 
9e-ph = SVA"^, which can be compared to the dimen- 
sionless conductance of the junctions. To cast this into a 
more understandable form, we write 



9e-ph/g - 



S V f A/kB 

'2-109 WK-V-3 2.8 • 10-21 m3 1^ 2.0 K 



12.9 kfi / ■ 



For example, for a small copper island and G = Gq, 
5e-ph/.9 = 1- When ge-ph/g <C 1 we can disregard the 
phonons altogether. 

For our purposes it is again enough to calculate 
Eq. (16) to the second order in ^. For phonons at zero 
temperature this can be done analytically. The electron- 
phonon heat current is given by 



rr _ ffe-ph rp5 

J^e-ph — 



and the heat current noise by 



Q _ ffe-ph rp6 

Finite temperature phonons add a term 



, ge-ph rp5 



(17) 



(18) 



(19) 



to the heat current. Due to the high powers of T these are 
negligibly small at low temperatures for not extremely 
large ge-ph/g- 



The average temperature as a function of bias volt- 
age, assuming zero-temperature leads and a vanishing 
electron-phonon coupling, is shown in Fig. 3. At low 
bias voltages the SINIS system exhibits overheating due 
to the second order Andreev processes. The sharp rise 
in temperature can be described analytically: Solving for 
T from ij = in the Hmit (A - V/2)/T > 1, F -I- 
yields the equilibrium temperature 



ln[(A/y)V^]- 



(20) 



As the voltage is increased, the cooling effect gets 
stronger, and the temperature reaches its minimum near 
V = 2A. We can again solve for the equilibrium temper- 
ature, which at y = 2A is given by 



T = 



7 



2/3 



r (1)53/2(0) 



(21) 



where 7 = V2TTrA/{gA) and 33/2(0) « 0.77. The mini- 
mum temperature is set by the magnitude of the Andreev 
processes. Above V = 2A the temperature is set mainly 
by the first order tunneling and is essentially independent 
of Ta- T(xV/2- A. 

We point out that instead of including the second order 
contributions the heating effect can also be reproduced 
by a finite broadening of the BCS density of states with 
the replacement e i-^- e + ir]'m the superconductor Green's 
function. The physical origin of this imaginary term can 
be, for example, the electromagnetic environment of the 
tunnel structure. 



B. Gaussian fluctuations 

The variance of the temperature, assuming a Gaussian 
distribution, can be estimated by^ 



Var(T) 



■'H 



2GthC 



(22) 



where C = '!t'^T/{36) is the heat capacity of the island 
and Gth is the heat conductance. The normalized vari- 
ance as a function of bias voltage in the Gaussian approx- 
imation is shown in Fig. 4. It is strongly peaked around 
the threshold voltage V « 2A, and falls rapidly when the 
voltage is tuned away from this threshold. 
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FIG. 3. (color online) Average temperature of the normal 
metallic island in a SINIS setup for some values of T . The 
inset shows that the scaled temperature t is independent of 
T near the scaled voltage « = 0. 



FIG. 4. (color online) Variance of the temperature of the 
normal metallic island in a SINIS setup for some values of T . 
The inset shows that the variance of the scaled temperature 
t is independent of T near the scaled voltage w = 0. 



Motivated by the 7^/^ dependence of the minimum 
temperature we define scaled variables i = T / {IS.')'^/^)^ 
V = {V/2 - A)/(A72/3), and d = 5/{Aj^/^), effectively 
eliminating 7 from the formulas in the limit 7^1. In 
terms of these scaled variables, we obtain 



Var(i) 



3d 



^5/2^3/2 (1^3/2(0) + agi/2{a) + a^g^i/2{a)) ' 

(23) 

where a = —v/t. For V — 2A, i.e., a = 0, we get using 
Eq. (21), 

Var(t) = ^. (24) 

Since the temperature is limited from below, the dis- 
tribution cannot be completely Gaussian. The low- 
temperature tail is always non-Gaussian. When the stan- 
dard deviation of the distribution approaches the expec- 
tation value, also the main body of the distribution be- 
comes non-Gaussian. Therefore, for d > (t) w 1, temper- 
ature fluctuations around V = 2A become non-Gaussian 
and of the same order of magnitude as the average tem- 
perature, rendering Eq. (23) incorrect. 



C. Induced current fluctuations 

Temperature fluctuations typically cause low- 
frequency fluctuations in the electric current. ^'^ The 
current noise due to temperature fluctuations is given 
by6 



Si. 



ind. 



/ dl/dT y 
KdH/dTj 



S 



(25) 



and is best characterized with the Fano factor of the SI- 
NIS system, Fsinis = Si/ {21). Near ^ = 2A we can 
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FIG. 5. (color online) Fano factor of the temperature fluctu- 
ation induced current noise for some values of J^. Inset shows 
the ^"^/"^ scaling of the Fano factor near the scaled voltage 
« = 0. 



use the analytical formulas for these quantities to obtain 
Si;ind./I ~ C7~^/^ (X where the numerical factor 
c « 0.17. This strongly peaked character of the Fano 
factor is shown in Fig. 5. For low the Fano factor can 
reach values of several hundreds, two orders of magnitude 
larger than the Fano factor -Fsinis ~ 1 due to intrinsic 
current fluctuations in the NIS junctions. 



D. Effect of non-zero temperature and 
electron— phonon coupling 

Non-zero temperature in the leads and/or an appre- 
ciable electron-phonon coupling have an effect on the 



6 




FIG. 6. (color online) Average temperature of the normal 
metallic island in a SINIS setup for T = 10~^, different values 
of Ts = Tph: 0.06A (blue), 0.08A (green), and 0.1 A (red), 
and different values of ge-p\i/g- (solid lines), 1 (dashed 
lines), and 10 (dash-dotted lines). The reference line for Ts = 
Tph = 0, Qe-pi^/g = is shown in black, which is lost under 
the solid blue line. 



expectation value and variance of the temperature on the 
island. The effect is strongest in the overheating regime, 
where the temperature of the island is highest. We are 
mostly interested in the effect around V = 2A, hovifever, 
since that is where the non-Gaussian effects should ap- 
pear. The effect to the expectation value of temperature 
in the scaled units around w = is shown in Fig. 6. Once 
the temperature of the leads approaches O.IA, the cool- 
ing effect begins to diminish. Similarly, once the phonon 
temperature approaches 0.1 A and the electron-phonon 
coupling parameter approaches unity, the cooling effect 
begins to diminish. The effect on the scaled variance is 
easily inferred from the effect on the expectation value: 
Since the variance is proportional to t~^/'^ , an increase of 
temperature by a factor of 2 around v — Q decreases the 
maximum variance by a factor of 0.35. 



IV. NON-GAUSSIAN STATISTICS 

A. Fokker— Planck equation 

To go beyond the Gaussian approximation we must 
turn our attention to the solution of Eq. (2). Using the 
equations derived in Subs. II B and II C the stationary 
state Fokker-Planck equation can be written in the form 







dE^H^^'^V{E)~dE 



-,-VdE 



-l)V{E) 



dE -ffe-ph7'(S) - dE 



(26) 



The terms on the first line originate from the tunneling 
of individual quasiparticles, second line describes two- 
electron Andreev processes, and the third line electron- 
phonon interaction. When the average energy on the 
island is large compared to voltage, we can expand e~^^^ 
to the second order in the Andreev term, and get 







dE { HV{E) - Oe 



S^ViE) 



(27) 



.(2) 



where H ^ H^i' + + H^-ph and , . 

•^ffe-ph- Using E — Tr'^T'^/{66) we can write a cor- 
responding equation for the probability distribution of 
temperature. This has the solution^ 



nt(T) = Cexp / dT 



\n^TH{T) - STOt (^#^) ' 



(28) 

where C is a normalization constant. In the limit Ta <C 



(1) 



A we can neglect S - 



He-ph, and 5*^ , 



-ph- 



These terms 



come with prefactors 7^^'^, 7^°/^, and 7^^^'^, respectively, 
whereas the other three terms have a common prefactor 
of 7. In the dimensionless variables introduced in the 
previous section the distribution simplifies to 



^5/2 



' 3d 



53/2(0) + 051/2(0) 



(29) 



For a = this becomes 
Vstit) = Ctexp 



2j.2 



~6d 



(30) 



which describes a distribution with an essentially Gaus- 
sian body and non-Gaussian tails. The condition of large 
average energy on the island translates to d ^ 1. For 
d > 1 Eq. (27) is no longer valid and we must take the 
higher derivatives into account in the Andreev term. 
Continuing in the limit Ta <SC A, and thus neglecting 



ph 



and 5* 



H,e-ph- 



we write 



h'-^^V{E) \ + Ta ( 



-.-Vdn 



l)V{E), 



= -dE \h^^^'>V{E)] + Fa [V{E -V)- V{E)] . (31) 

Its solution depends on V and Ta, resulting in six qual- 
itatively different regimes that are studied in detail in 
Ref. 8. In general, for F^ < gA^^^S'^^'^, or, equivalcntly, 
J" < (5/A)3/4, and IF - 2A| < VSA, the resuhing prob- 
ability distribution is non-Gaussian. 



B. Langevin equation 

The Fokker-Planck equation, Eq. (31), can be equiv- 
alcntly presented in the form of a Langevin equation. 
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Starting from the partition function, Eq. (1), and using 
S{^,E) = + rA(e«^ - 1) as the action, we can 

expand the exponentiated action in powers of F^: 



^^^^ E ^ n / d^^e- 



m=0 k=l' 



X exp ■ 



dt£_ 



(32) 



By performing now the functional integral over ^ we 
end up with a delta functional of the quantity in square 
brackets. This means that only paths which satisfy the 
Langevin equation 



(33) 



contribute to the partition function. The energy is in- 
creased in increments of V due to Andreev events at 
random times tk, and between the Andreev events en- 
ergy changes deterministically as dictated by the energy 
flow ij^^-* due to single particle processes. The Andreev 
events are Poisson distributed as can be seen from the 
statistical weighing factor in the functional integral. 

Using the Langevin equation we can simulate timelines 
of energy, which are sometimes more informative than 
the full distributions obtained from the Fokker-Planck 
equation. These timelines are presented in Ref. 8. 



C. Non-Gaussian current fluctuations 

The electric current through the SINIS structure near 
y = 2A is given by^ 



(34) 



where e = E/{2A), Tt = ViMA/tt, u = {V/2 - A)/Tt, 
and = gT^^^ /{tt^/2A). A fluctuating internal en- 
ergy (or temperature) leads directly to a fluctuating 
electric current, and the probability distribution for the 
fluctuations of the electric current can be found by a 
change of variables from the probability distribution of 
energy. Three distributions in the strongly non-Gaussian 
regime are shown in Fig. 7 for the dimensionless current 
J = T^Ttl / A. The current distributions mirror the non- 
Gaussian features of the distributions for internal energy 
on the island. For larger u, current increases exponen- 
tially, and the distribution turns into a Gaussian peak 
at larger j. For smaller u, the distribution is a narrow 
Gaussian peak near j = 0. 




FIG. 7. (color online) Probability distribution for the dimen- 
sionless current j = TrTtI / A in the non-Gaussian regime, 
r,4T"r — 0.5, for three different values of bias voltage. The 
sharp peaks in the two leftmost distributions correspond to 
the case where the internal energy on the island is E — 2A. 



V. DETECTION OF TEMPERATURE 
FLUCTUATIONS 

A possible way to detect the temperature fluctuations 
is by detecting the induced fluctuations in some other ob- 
servable. As shown in Ref. 8, the statistics of temperature 
fluctuations can be gathered by monitoring the current 
through the device in real time. The proposed measure- 
ment scheme is shown in Fig. 8. The temperature de- 
pendent current through the SINIS structure, Iswis{T), 
generates a fluctuating voltage Vx across a shunt resis- 
tor i?shj which is then amplified and detected in time 
domain, with a fast oscilloscope, for example. 

In order for the temperature fluctuation induced noise 
to be detectable, the additional noise generated by the 
shunt resistor and the amplifier must be smaller than 
the actual signal to be measured. This implies that 
■\/ variVx) <C (Vx), where the overbar denotes a time av- 
erage [see Eq. (35) below]. To make further progress we 
decompose the fluctuating voltage into a signal part and 
a noise part, Vx = V + SV, and average over a short mea- 
surement period r ^ t^, being the relaxation time for 
the temperature fluctuations: 



I 



dt{Vx{t))^V. 



(35) 



Variance of the measured signal is given by 
{Vx')-{Vxr=^ r f dm'{5V{t)5V{t')), 



1 



Jo 

T pT 



dtdt' 



dio 
2.' 



(36) 
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10 kn, Rt = 4.3 kn, a = 2 K X feg, .5 = 
and r = 10 mK, we get {V^) / y/vai{V^) 



10-4 K X fcs, 
= 8, implying 



that the temperature fluctuation induced noise should be 
readily detectable. 



z 



\ 
/ 



T 



Rsh 



FIG. 8. Example measurement scheme for the detection of 
temperature fluctuations via current fluctuations. The cur- 
rent through the SINIS structure at a fluctuating tempera- 
ture T, /siNis(T), generates voltage fluctuations across the 
shunt resistance Rsh- These are amplifled and recorded in 
time domain. 



where Sy^ (w) is the spectral noise power. Typically, the 
dominating source of extra noise is the thermal noise in 
the shunt resistor. In this case Sy^ is independent of 
frequency^''^^ and the variance becomes 



var(V;) = 



dyn 



T(i?sh + i?dyn)^ 



(37) 



where -Rdyn = (Vbias — V)/I is the dynamical resistance 
of the SINIS structure. 

The signal-to-noise ratio becomes 



Vvar(y^) 




1/2 



(38) 



whore on the second line we have used r = krr and the ex- 
pressions from Ref. 8. Using the numbers = 0.1, -Rsh = 



VI. CONCLUSIONS 

Building on the recently developed theory of tempera- 
ture fluctuation statistics, we have studied the tempera- 
ture fluctuations and the associated current fluctuations 
in a SINIS structure. We have shown how the Andreev 
processes lead to a heating of the normal metal island 
and how the cooling effect near V = 2A is accompa- 
nied by fluctuations in the island temperature. We have 
also shown that these temperature fluctuations lead to 
current fluctuations with a Fano factor that is large com- 
pared to the Fano factor of the intrinsic shot noise. 

In contrast to the case of the overheated single-electron 
transistor,^ the observation of the temperature fluctua- 
tions docs not require tunnel junctions of very high re- 
sistance. The statistics of temperature fluctuations can 
be gathered by monitoring the electric current through 
the device in time domain, and extracting the statistics of 
temperature from the statistics of the fluctuating electric 
current. 

In addition, the cooling effect of the SINIS structure 

ensiires that the island temperature near V = 2A is low 
enough that the electron-phonon coupling does not mask 
the fluctuations. Due to the poor thermal conductivity 
of superconductors, the practical challenge is rather to 
keep the superconducting reservoirs in equilibrium and 
at a low temperature. 
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Appendix A: Derivation of the electron-phonon action 



The electron phonon action is equal to the sum of all one particle irreducible electron-phonon vacuum diagrams. 
The lowest order diagrams are shown in Fig. 9. It is enough to consider the diagram on the left, tlie diagram on the 
right vanishes due to zero momentum exchange. In terms of Keldysh Green's functions we have 

5e-ph = Tr/^^yE^(k, e)7|5k-kf £>(k - k', e - e')G{k' , (Al) 

k.k' 
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k',e' 

FIG. 9. Two lowest order electron-phonon diagrams contributing to the effective action. 



where Tr is trace over Keldysh and Nambu indices. The absorption and emission vertices are, respectively/® 

lio=h^ok, lii=hrfl (A2) 

Since the electrons participating in the interactions lie close to the Fermi surface, we can use the quasiclassical 
approximation to write 



J d^^d^FG(ek,e)7|5fcWe.-e,,)l'^(fci='(ek-ek'),e-e')G(ek',e')7, 



(A3) 



where dSp is an area element of the Fermi surface, u the density of states at the Fermi surface, and ek unit vector in 
the direction of k. 

Assuming that phonons are not affected by the electrons, i.e., that they are in equilibrium at temperature Tph, 
= {D^ — D^) coth(e/2Tph). For the electron Green's function we use Eq. (3). Carrying out the matrix and tensor 
operations results in 

X [e«(-^')„(e - e') (/(e')(l - /(e)) + (1 - /(-e'))/(-e)) 
+e-«(-^') (1 + n(e - e')) (/(e)(l - fie')) + (1 - fi-e))f{-e'))_ 

- e«(-^')n(e - e') [{/%){! - /(e)) + /^(e)/(-6)) (.f«(e')(l - /(-^')) + /^(e')/(e')) 
+ if''{e)f{-e) + /■^(e)(l - /(e))) {.f%').fie') + /^(e')(l - f{-e')))] 

_ e-«-^') (1 + n{e - e')) [{f%).fie) + /^(e)(l - /(-e))) {.f^{e').f{-e') + /^(e')(l - /(e'))) 

+ (/«(e)(l - /(-e)) + /^(e)/(e)) {f%'){l - f{e')) + /^(e')/(-e'))] } (A4) 

where 



If /(— e) = 1 — /(e) this simplifies to 

de de' Vv^ 



'^-P^^/ / <iSFdS'^\g,,^,,_,^,)f{D^ - D^){N{e)N{e') - F{e)F{e')) 

X (e«-^')„(e - e')/(e')(l " /(e)) + e'^^-^'^ (1 + n(e - e')) /(e)(l - /(e'))) , (A5) 



where 

Asgn(e) 



F(e) = Re 



- A2y' ■ 

Let us now concentrate on the spectral function of phonons. In equilibrium 

A(q, e) = i (i?«(q, e) - i?^(q, e)) = 27r ((5(e - £q) - 6{e + s^)) , (A6) 
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where Sq = wr>|q| is the (acoustical) phonon dispersion relation with the speed of sound vd- In the jellium model the 
e-ph coupling is j^qp = where g is now some dimensionless coupling constant. For a spherical Fermi surface, 

and any function /, 

J dSpdS'pfileu - ek' I) = S^r^ ^ dkkf{k). (A7) 

Using this result we get 



/ 



d5fd5^.|5fc^(e,_e^,)|'^(fcF(ek-ekO,e-e') = ^^{e - e'f sgn{e - e')9{2e^,^ - \e - e'\), (A8) 



"max 



where £max = vokp- This can now be inserted to Eq. (A5), which after identifying S = 2A-Kg'^v'^C,{b) / e^^^^ leads to 
Eq. (16) for a normal metal island. 
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